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Abstract 



The construction and use of a dual radiator Ring Imaging Cerenkov (RICH) detec- 
tor is described. This instrument was developed for the HERMES experiment at 
DESY which emphasises measurements of semi-inclusive deep-inelastic scattering. 
It provides particle identification for pions, kaons, and protons in the momentum 
range from 2 to 15 GeV, which is essential to these studies. The instrument uses 
two radiators, C4F10, a heavy fluorocarbon gas, and a wall of silica aerogel tiles. 
The use of aerogel in a RICH detector has only recently become possible with the 
development of clear, large, homogeneous and hydrophobic aerogel. A lightweight 
mirror was constructed using a newly perfected technique to make resin-coated 
carbon-fiber surfaces of optical quality. The photon detector consists of 1934 photo- 
multiplier tubes (PMT) for each detector half, held in a soft steel matrix to provide 
shielding against the residual field of the main spectrometer magnet. 
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1 Introduction 

The HERMES experiment [1] is a study of 
the spin structure of the nucleon which em- 
phasizes on an unambiguous measurement 
of pion, kaon, and proton semi-inclusive 
spin asymmetries in deep-inelastic scatter- 
ing (DIS). These asymmetries provide in- 
formation on the flavor dependence of po- 
larized structure functions and the sea po- 
larization. However, most of the hadrons 
produced in HERMES [2] lie between 2 and 
10 GeV, a region in which it has not been 
feasible to separate pions, kaons, and pro- 
tons with standard particle identification 
(PID) techniques. Ring imaging Cerenkov 
(RICH) and threshold Cerenkov systems 
using heavy gases [3], such as C 4 F 10 , at at- 
mospheric pressure are useful only for en- 
ergies above 10 GeV since the kaon 
threshold for Cerenkov radiation is typi- 
cally higher than 9 GeV. Because of sub- 
stantial multiple scattering and brems- 
strahlung, the use of a high pressure gas 
system is not technically feasible in HER- 
MES. Clear liquid radiators are only use- 
ful for hadron identification below roughly 

2 GeV because of their very low Cerenkov 
light thresholds and large chromatic dis- 
persion. 

With the recent development [4,5] of new 
clear silica aerogel with a low index of re- 
fraction, this difficult energy region can now 
be spanned. Aerogels have long been used 
in threshold Cerenkov counters [6] , but their 
use as a radiator in a RICH system is a 
recent development [7] which has resulted 
from the availability of material with ex- 
cellent optical properties. This paper re- 
ports the successful use, for the first time, 
of clear aerogel in combination with a heavy 
gas, C4F10, in a RICH detector. This dual- 
radiator RICH detector provides clean sep- 



aration of pions, kaons, and protons over 
most of the kinematic acceptance of the 
HERMES experiment. Such a configura- 
tion was first proposed for the planned 
LHCb experiment [8]. 

The following section 2 describes the gen- 
eral design of the HERMES RICH detec- 
tor and the properties of its key compo- 
nents. Section 3 describes briefly the on- 
line monitoring of the system, section 4 the 
reconstruction and particle identification 
and section 5 the necessary software cor- 
rection of the mirror alignment. Section 6 
finally is dedicated to the performance of 
the detector. 



2 Detector Design 

2.1 Design Requirements 

HERMES is located in the East hall of 
the HERA storage ring complex. The spec- 
trometer is split into two halves above and 
below a central horizontal magnetic shield- 
ing plate through which the HERA beams 
pass [1]. As a result, the RICH detector 
also consists of two symmetric RICH mod- 
ules, the top and the bottom RICH sec- 
tions. The RICH is positioned between the 
two rear tracking chambers, so that the 
free space is only 1.25 m deep. The outer 
body and the gas control system were taken 
unchanged from the previous threshold 
Cerenkov counters. 

Figure 1 shows the hadron momentum spec- 
tra in the HERMES acceptance obtained 
from a GEANT based Monte Carlo (MC) 
simulation of the deep inelastic scattering 
of 27.5 GeV positrons off a hydrogen tar- 
get. The fall-off of the spectra at low mo- 
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Fig. 1. Monte Carlo hadron momentum spec- 
tra within the HERMES acceptance. 

menta is due to the field of the spectrom- 
eter magnet, which severely limits the ac- 
ceptance at lower momenta. About 95% of 
all hadrons in the acceptance are found in 
the range of 2.0 to 15.0 GeV. This defines 
the momentum range over which clear par- 
ticle identification should be provided. 

The low end of this range determines the 
index of refraction necessary for the aero- 
gel. A value of n(A=633 nm)=1.03 was cho- 
sen since it leads to a kaon threshold of 
2 GeV. The Cerenkov angles produced by 
the combination of this aerogel and the 
heavy gas (C4F10) for pions, kaons and pro- 
tons are plotted in figure 2 as a function 
of particle momentum. The corresponding 
threshold momenta are listed in table 1. 
All pion momenta within the spectrometer 
acceptance are above the pion threshold 
momentum for aerogel of 0.6 GeV, 90% of 
the kaon and 78% of the proton momenta 
are above the kaon threshold of 2.0 GeV. 

The high end of the momentum range fixes 
the number of photons that must be de- 
tected for full hadron separation. The pa- 



Fig. 2. The Cerenkov angle 9 versus hadron 
momentum for the aerogel and C4F10 gas ra- 
diators. 





aerogel 


C4F10 


n 


1.0304 


1.00137 




4.03 


19.10 


7T 


0.6 GeV 


2.7 GeV 


K 


2.0 GeV 


9.4 GeV 


P 


3.8 GeV 


17.9 GeV 



Table 1 

Cerenkov light thresholds for pions, kaons and 
protons. The index of refraction n is given at 
633 nm, fit = 1/n is the threshold velocity 
and 7 t = 1/ y/l 



Pi 



rameter to be considered is p ma x , the max- 
imum separation momentum [9]. This is 
defined as the maximum momentum for 
which the average photon emission angle 
of two particle types (with masses mi and 
1712) is separated by a number of standard 
deviations n„ : 
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where kf — tan 6 ■ ctq / y N is the RICH de- 
tector constant, N is the number of sepa- 
rately detected photons, 6 is the Cerenkov 
angle and ag the standard deviation of the 
reconstructed photon angle distribution. In 
the design, n a = 4.652 was chosen, as it 
corresponds to a misidentification of the 
particle in 1% of the cases, assuming equal 
fluxes for the two particle types, an average 
detector response (in yield and resolution) 
and no background. 

Assuming <jq to be 7 mrad (see tables 3 
and 4) it follows from (1) that p max {ji,K) 
= 15 GeV requires N for the gas to be 
12. This requirement leads to the design 
values for p ma x given in table 2. In this es- 
timate it was assumed that the number of 
separately detected photons from the aero- 
gel is 10. Figure 3 illustrates the overlap 
between the momentum regions for both 
radiators. The lightly shaded region indi- 
cates where the particle can be identified 
based on whether or not a ring is present 





aerogel 


C4F10 


*/ 

Pmaxi^/K) 
Pmax(K/p) 


5.46 • 10" 4 
6.7 GeV 
11.2 GeV 


1.07 • 10~ 4 
15.0 GeV 
25.3 GeV 



Table 2 

Maximum separation momenta p n 




effective threshold 



Fig. 3. Momentum ranges for hadron sep- 
aration in aerogel and C4F10. Between the 
dashed lines the hadrons can be separated. 



at all. In this region the detector acts like a 
threshold Cerenkov. In the darkly shaded 
region the identification is based on the av- 
erage reconstructed angles. The plot con- 
siders each radiator separately, but the PID 
algorithms will combine the information 
from the two. The momentum region for 
which the identification of pions, kaons and 
protons is possible is limited by the kaon 
threshold momentum for aerogel at 2.0 GeV 
and by the maximum separation momen- 
tum for 7r/K separation in C4F10 at 15.0 
GeV. 



2.2 General Design Parameters 



The geometry which was adopted for the 
Cerenkov radiators and ring imaging sys- 
tems is shown in figure 4 [10]. The body 
of the counter is constructed of aluminum, 
with entrance and exit windows made of 1 
mm thick aluminum. The volume of each 
half is approximately 4000 1. The size of 
the entrance window is 187.7 cm by 46.4 cm 
and the exit window 257.0 cm by 59.0 cm. 



PMT plane 



Dry N : 



Particle 




Aerogel Tiles 



Fig. 4. Basic geometry and radiator con- 
figuration for the HERMES dual radiator 
RICH (not to scale). 
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aluminum box 



Fig. 5. A cutaway schematic view of the (top) RICH counter. 



A gas control system recirculates the radi- 
ator gas through the main volume, keep- 
ing the gas at a slight overpressure with 
respect to atmosphere. The aerogel radia- 
tor is an assembly of tiles configured to fill 
the entrance of the detector with an aero- 
gel thickness of 5.5 cm. The unoccupied 
volume of the detector behind the aerogel 
is filled with the gas radiator, C 4 Fi . A 
spherical mirror array located at the rear 
of the radiator box images the Cerenkov 
light cones onto a focal surface located 
above (below) the active volume. 

The radius of curvature of the mirror ar- 
ray is 2.20 m. It was chosen to give a fo- 
cal surface location in the accessible re- 
gion above (below) the forward region of 
the radiator boxes and to provide a de- 
tector plane of tractable dimensions. The 
optical axis of the array, the perpendicu- 
lar to the mirror surface at the center of 
the array, is inclined at an angle of 26 de- 
grees to the horizon. The photon detector 
is located outside of the mirror optical axis 
with its axis inclined at an angle of 40 de- 
grees to the horizon so as to intercept the 
mirror surface at a distance of 90 cm. The 



focal length of the mirror is 110 cm. The 
boxes are fitted with gas connections and 
pressure regulators which provide a con- 
tinuous controlled flow of recirculating gas. 
An open section of one of the RICH coun- 
ters is shown in figure 5. 

The size of a useful detector surface was 
evaluated by an MC simulation which in- 
cluded an early version of the RICH ge- 
ometry described above. The simulation 
showed that 95% of the centers of the rings 
and 90% of all the photons are contained 
in a planar surface 60 cm high and 120 cm 
wide (0.72 m 2 surface area). These dimen- 
sions were used as lower limits in the final 
design of the photon detector. 

The inner walls of the box are blackened to 
reduce wall reflections. An array of green 
light-emitting diodes (LED) is installed to 
provide test and calibration pulses for the 
photon detector. They are located on the 
face of the mirror, so as to provide an ap- 
proximately uniform illumination of the 
photon detector surface. 

As explained below, most of the useful pho- 
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tons contained in rings emitted from the 
aerogel are in the visible light region. A 
good choice for the photon detector is then 
an array of photomultiplier tubes (PMTs). 
In the present design 3/4 inch tubes were 
used, leading to a pixel size of 23.3 mm. A 
total of 1934 PMTs is necessary to cover 
the photon detector area. 

Table 3 lists estimates for the different con- 
tributions to the single photon resolution. 
A9 em results from the uncertainty in the 
emission vertex along the track in the radi- 
ator. The detector granularity determines 
A9 pix , which follows from the size of the 
"pixels". These two contributions affect 
both the aerogel and the gas angles. The 
chromatic aberration A6 c h r derives from 
the variation of the index of refraction with 
respect to the wavelength and is only im- 
portant for the resolution in the aerogel. 
These quantities do not add quadratically 
and are only listed separately to indicate 
their relative importance. Their combined 
effect is calculated in the MC and is listed 
as A9 M c- The resolution is dominated by 
the pixel size, which is determined by the 
size of the PMT and how much material 
must surround it to provide adequate mag- 
netic shielding. 





aerogel 


C4F10 


A6> em 


1.8 mrad 


2.2 mrad 


AQpix 


5.6 mrad 


5.2 mrad 


A0 chr 


2.5 mrad 






7.1 mrad 


7.2 mrad 



Table 3 

Individual contributions to the single photon 
resolution in aerogel and C4F10. 



2.3 Aerogel Radiator 

New production techniques [4,5] have yiel- 
ded aerogels with much more uniform and 
smaller colloidal SiC>2 particle structure. 
The more uniform granular structure sup- 
presses Rayleigh scattering of visible light 
and greatly enhances its transmission. Ma- 
terial with a refractive index in the range 
1.01 — 1.10, suitable for use in the mo- 
mentum range of interest in HERMES, is 
now commercially available from Matsu- 
shita Electric Works 1 in the form of tiles 
of average dimensions 114 by 114 by 11.3 
mm 3 . The optical properties of the aero- 
gel tiles have been studied in detail using a 
sample of 200 tiles. The transmittance was 
measured for each tile for the wavelength 
range of 200-900 nm. Typical values for in- 
dividual tiles at 400 nm are around 0.67. 
The transmittance T was parameterized 
as a function of wavelength A and aerogel 
thickness t by the Hunt formula [13]: 

T = A-e^ (2) 

The measurements can be described well 
with a value of C ■ t — 0.0094 /im 4 and 
A = 0.964. More details on the optical 
properties of the aerogel tiles can be found 
in [11-14]. 

Out of 1680 tiles produced in a single pro- 
duction batch, 1020 tiles were judged to 
be usable for the HERMES RICH. The se- 
lection was based on the measured refrac- 
tive index and the sizes of the individual 
tiles [11]. Of the 1020 tiles, 850 were ac- 
tually used. For those tiles, the average in- 
dex of refraction is 1.0304, with a spread of 

1 Advanced Technology Research Labora- 
tory, 1048 Kadoma, Kadoma-shi, Osaka-fu, 
Japan 571, contact: M. Yokoyama, email: 
yokoyama® crl. mew .co.jp 
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3.6 • 1CT 4 [12,11]. The tiles are stacked in a 
container consisting of an aluminum frame 
with a 1 mm aluminum entrance window 
and a 3.2 mm UVT-lucite exit window. 
To guard against possible degradation of 
the aerogel by the C4F10 environment, the 
aerogel container is sealed gas tight and 
dry nitrogen is continuously circulated 
through the box at a slow rate. 

The aerogel tiles are stacked in 5 layers, 
with 5 horizontal rows, and 17 vertical co- 
lumns as required to span the spectrom- 
eter acceptance. Black plastic spacers of 
appropriate thicknesses between the alu- 
minum frame and the tiles prevent them 
from shifting while the radiator is moved. 
The aerogel wall thickness was chosen as 
the optimal point for maximizing the un- 
scattered light yield relative to scattered 
background photons [13,14]. The tiles are 
stacked according to their measured refrac- 
tive index, thickness, and surface quality. 
To achieve the best ring resolution, tiles 
with similar refractive indices are placed 
together in the same stack, so that par- 
ticles passing through the radiator emit 
photons with very similar Cerenkov angles. 
Opaque black sheets of tedlar between the 
aerogel stacks reduce distortions by absorb- 
ing photons that cross stack boundaries. 

The choice of the material for the exit win- 
dow was driven by its transmission proper- 
ties. Due to the proportionality to A~ 4 of 
the Rayleigh scattering cross section, the 
low wavelength end of the light yield spec- 
trum is dominated by rescattered photons. 
Since they only contribute to the back- 
ground, the exit window material was se- 
lected to absorb most of them. A good 
match was found in ultraviolet-transmitting 
(UVT) lucite, which has an absorption cut- 
off of 290 nm [14] (see figure 6). The thick- 
ness of the UVT window, 3.2 mm, was 



chosen to be as small as possible while 
still providing sufficient mechanical stabil- 
ity. It contributes less than 1% of a radia- 
tion length. The Cerenkov light generated 
in the lucite itself is emitted at such high 
angles that it is effectively trapped by total 
internal reflection for most forward tracks. 



2.4 Mirror Array 

The mirror arrays each consist of eight seg- 
ments that are mounted in two rows of four 
(see figure 5). Both the mirrors and the 
mounting frame were fabricated by Com- 
posite Mirror Applications (Tucson, Ari- 
zona, USA) 2 . The backing of the segments 
is fabricated from a graphite fiber compos- 
ite [15] coated with an epoxy film yield- 
ing an optically smooth surface. Surface 
roughness was specified to be less than 5 nm 
The surface is aluminized to provide a re- 
flectivity above 85% for light in the 300- 
600 nm range (see figure 6). The segments 
are held on a carbon-fiber frame with indi- 
vidual adjustable three point mounts. The 
dimensions of each array are 252.4 cm by 
79.4 cm; each weighs less than 13 kg. This 
guarantees that the mirrors contribute neg- 
ligibly (less than 1.0% of a radiation length) 
to the total radiation length of the detec- 
tor, which is dominated by the freon gas 
(5%), the aluminum windows (3%) and the 
aerogel (2.8%). 

The short time scale of the RICH installa- 
tion (13 months from approval to installa- 
tion) only allowed for a simple alignment 
procedure when the mirrors were installed. 
A point light source at the nominal radius 
of curvature of the array was imaged on a 



2 web: http://home.earthlink.net/~bromeo/ 
contact: R. Romeo, email: bromeo@earthlink.net 
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Fig. 6. Efficiencies for transmission and detection of aerogel (left) and gas photons (right). The 
included contributions are the mirror reflectivity (mirror), the PMT cone efficiency (cone), the 
effect of the air gap between PMT and quartz window (loss), the light transmission through 
lucite and the quantum efficiency of the PMT bialkali cathode. The sharp edge in the cone 
efficiency results from the sharp cutoff in reflectivity at 200 nm. The shaded histogram is the 
total efficiency to detect an emitted photon. 



screen. The mirror segments were then ad- 
justed to minimize the size of the reflected 
image. To neutralize the unexpected flex- 
ibility of the mirrors, their mounting was 
fine-tuned by adjusting additional pressure 
points where necessary. A nominal 90% of 
the reflected light was contained in a spot 
corresponding in size to the pixel size of 
23 mm. The actual mirror alignment was 
determined from the data itself (see sec- 
tion 5). 



2.5 Photon Detector 



As discussed, the useful photons emitted 
in the aerogel are largely limited to visible 
wavelengths. The radiation emitted in the 
gas does have a component extending into 
the UV region, but since the photon yield 
from the gas is sufficiently high, the detec- 
tor response in the UV is not crucial. Con- 
sequently, photon-to-electron conversion by 
photocathodes typical of commercial pho- 
tomultipliers provides a simple and robust 



photon detector technology. 

While photomultipliers with a diameter of 
0.5 inch are commercially available, their 
cost and the time of construction for two 
large 0.5 inch arrays were unacceptable. 
Instead, Philips XP1911 photomultipliers 
were chosen, with a diameter of 18.6 mm 
(0.75 inch) and a guaranteed minimum ac- 
tive photocathode diameter of 15 mm [16]. 
This led to a pixel size of 23.3 mm. A9 pix 
can be estimated at the center of the de- 
tector by assuming that both the radia- 
tor and the detector are at the focal plane. 
This gives A9 P i X = s/ (4/) = 5.3 mrad (for 
circular pixels), where / is the focal length 
(110 cm) and s is the linear pixel size. How- 
ever, this is a simplifying assumption and 
an MC simulation leads to the values in 
table 3. 

The XP1911/UV (green enhanced) PMT 
was chosen for its broad quantum efficiency 
curve that matches the Cerenkov light spec- 
trum of the aerogel well. Since the effi- 
ciency extends into the UV, it is also a 
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good choice for the gas photons, as is clear 
from figure 6. This figure shows the differ- 
ent contributions to the total photon de- 
tection efficiency for both radiators. 

An automated PMT testing procedure was 
used to select PMTs with a broad plateau. 
High voltage (HV) plateau curves were mea- 
sured for each tube with a 470 nm light 
emitting diode pulser. In addition, a low 
dark current single photoelectron pulse rate 
was required. After the installation of the 
RICH in the HERMES spectrometer a noise 
rate of about 1 kHz was reached, corre- 
sponding to 1 PMT firing per detector half 
every five events (gate length 100 ns). The 
PMTs were sorted according to their gain 
characteristics and their noise level. With 
an especially developed low current (0.040 
mA at 1300 V) high voltage divider it was 
possible to supply 32 PMTs from a single 
HV channel. The divider also limited the 
heat generated in the base of the PMTs, 
thereby significantly reducing the cooling 
requirements of the detector. 



184.15. 
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slip-fit insert 
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Fig. 7. Schematic photon detector design. All 
units are in mm. 

vidual PMT cavities was provided by thin 
quartz windows which were glued to soft 
steel cone-shaped inserts that in turn were 
glued into the inner face of the PMT ma- 
trix plate. 



The photomultipliers were arranged in a Only 38% of the area of the focal plane 
hexagonal closed packed matrix of 147.4 cm is covered by the PMT photocathodes. An 
by 62.8 cm, centered on the mirror focal aluminized plastic foil funnel was inserted 
point. This configuration was first used by into each soft steel entrance cone to in- 
the experiment E781-SELEX at FNAL [17]. crease the coverage to 91% of the pho- 
The elementary cell of the array is a hexagon ton detector surface and to improve the 
with the photomultiplier at the center. A reflectivity of the cone surface. The fun- 
light-gathering cone is used to increase the nels extend beyond the steel plate with 
photon collection efficiency by reducing the an opening diameter of 23.3 mm, so mini- 
dead space between the photocathodes. The mizing the dead space between the PMTs. 
distance between two adjacent cell centers They provide a high reflectivity above A = 
is 23.3 mm, yielding an elementary cell sur- 200 nm. Their effect on the geometric col- 
face of 426 mm 2 . The 1934 PMTs in each lection efficiency is shown in figure 8 as a 
detector half are arranged in 73 columns function of the angle between the incom- 
of alternately 26 and 27 PMTs each. A ing photon and the normal to the detec- 
section of the soft steel matrix plate and tor plane (incident angle). The influence 
PMT packing design is shown in figure 7. of the funnels is found from the difference 
The matrix plate provides the gas seal for between the total efficiency and the effi- 
the C 4 F 10 volume. The gas seal for indi- ciency for "direct" detection. 
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Fig. 8. Collection efficiency of the funnels ver- 
sus incident angle (solid points) from a Monte 
Carlo simulation. The other curves show the 
components of the efficiency in terms of the 
number of "bounces" , i.e. reflections from the 
foil surface. 

The stray fields from the HERMES spec- 
trometer magnet at the location of the pho- 
ton detector matrix can reach up to 90 G 
perpendicular to the detector plane, requir- 
ing careful magnetic shielding of the PMTs. 
Individual PMTs were wrapped in a 100/im 
thick /x-metal sheet and the matrix plates 
holding the PMTs were constructed from 
high permeability (C-1008) steel. The soft 
steel matrix, combined with the //-metal 
shielding and the soft steel inserts, reduces 
the magnetic fields to negligible levels, which 
guarantees that the PMT gains are not sig- 
nificantly affected. 



2. 6 Detector Readout and Event Format 

The readout of the photon detector is per- 
formed by the commercial LeCroy PCOS4 
acquisition system, upgraded for the HER- 
MES application. Each detector half is read 



out by a set of 8 symmetric backplane sec- 
tions, each housing 15 or 16 PCOS4 cards. 
Each card in turn processes signals from 16 
PMTs. This reduces the number of cables 
needed to 10 twisted pairs cables. Only 
digital information - when the pulse ex- 
ceeds the threshold of 0.1 photoelectrons 
- is recorded. The system is characterized 
by high input sensitivity (the threshold is 
3000 electrons) and high amplification 
(—4.3 /iV per electron). Because of the 650 
W power dissipated per detector half, each 
RICH section has a cooling system inte- 
grated into the RF-shield and electronic 
enclosure. A closed stream of cooling air 
(360 m 3 /h) is guided over the PCOS4 elec- 
tronics on both sides of the PM-matrix 
and over the PMT matrix itself. The heat 
transferred to the air flow is removed by a 
water-cooled heat exchanger. 

The data from both RICH detector halves 
are treated in the same manner as spec- 
trometer chamber data. The PCOS4 sys- 
tem generates a data stream on an event- 
by-event basis. The RICH data consists of 
a RICH hit table for each detector half, 
where a PMT 'hit' simply refers to a PMT 
that fired during the event time window. A 
RICH mapping file is used to link the data 
channel numbers computed in the decod- 
ing process to the PMT location in the de- 
tector matrix. This information can then 
be used to generate the spatial hit coor- 
dinates and the hit pattern in the focal 
plane. 



3 Online Monitoring 

HERMES uses a client-server based sys- 
tem for the monitoring of detector operat- 
ing parameters [18,19]. In the case of the 
RICH this includes high and low voltages, 
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Fig. 9. Two-dimensional tubemap of relative 
yield for the bottom RICH detector. 

gas pressure and composition, and the tem- 
peratures at the photon detectors. The in- 
formation is available in the form of Tcl/Tk 
based displays [20]. 

In analogy to wiremaps for the tracking 
chambers, "tubemaps" are used to moni- 
tor the performance of the PMTs, the high 
voltage and the electronics. An example of 
the standard two-dimensional overview is 
shown in figure 9. The increased hit den- 
sity in the center is due to the gas rings. 
This plot shows that more than 99% of the 
tubes are functional. More detailed maps 
are available for diagnosing operational prob- 
lems. 

The RICH data have also been integrated 
into the online HERMES event display sys- 
tem to allow the online analysis of events. 
Figure 10 shows the standard view for the 
case of a clean single electron event. The 
event display combines top and side views 
of the spectrometer with views of the top 
and bottom RICH photon detectors. It rep- 
resents the grid of 1934 PMTs in the form 
of dots. All PMTs that fire are marked by 
filled circles. In general particles with mo- 
menta above both the aerogel and C 4 F 10 
thresholds will generate roughly concentric 
rings: a small inner ring resulting from the 
photons generated in the gas and a larger 
outer ring formed by the aerogel photons. 



Electrons and positrons that are detected 
in the RICH always have momenta above 
both radiator thresholds. They should there- 
fore always have both rings present, which 
makes them ideal for the online monitor- 
ing of the data quality. The average num- 
ber of PMT hits for single lepton events in 
a run is used as the monitoring parameter; 
it is usually about 20. Figure 11 shows this 
parameter in a typical online data quality 
plot. 
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Fig. 10. Online display of a single event. 
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Fig. 11. Average number of hit PMTs per sin- 
gle lepton event plotted versus run number. 



4 Reconstruction and Particle Iden- 
tification 



4-1 Principles of the Reconstruction 



Since the sensitive face of the flat photon 
detector does not conform to the true mir- 
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ror focal surface, the detected 'rings' are 
not circles, but in general asymmetrically 
distorted ellipses. Hence a ring fitting pro- 
cedure is not feasible for an accurate re- 
construction of the Cerenkov angles. The 
raw data that provide the input for any re- 
construction procedure are the list of PMT 
hits per event and the parameters of the 
particle tracks determined by the tracking 
detectors [1]. 

In the case of the aerogel rings on average 
only about 10% of the PMTs that lie on 
the aerogel ring will be hit. Hence mul- 
tiple hits are infrequent and each PMT 
hit represents one photoelectron. The av- 
erage number of hits on the aerogel ring of 
an ultra- relativistic (/3 ~ 1) particle that 
suffers no acceptance effects (aerogel tiles, 
mirror) is about 10. The sparsity of the 
hits on the aerogel ring makes the detec- 
tion of this ring susceptible to background, 
noise and acceptance problems. The situa- 
tion for gas rings is very different. The av- 
erage multiplicity of photoelectrons per hit 
for the gas rings is about 2 (based on MC 
simulations). This might lead to the as- 
sumption that all PMTs that lie on the gas 
ring will be hit. However, due to Poisson- 
statistics this is not correct. On average 
about 23 PMTs are touched by the gas ring 
of an ultra-relativistic particle, but only 12 
of those PMTs actually fire. Most of them 
are hit by several photo-electrons. 

Figure 12 shows the offline RICH event 
display for a three-track event that illus- 
trates several typical features of HERMES 
RICH events. The event shown has a 14.6 
GeV electron 3 and a 1.5 GeV pion in the 



3 Electrons and positrons are identified by 
a combination of calorimeter, preshower and 
TRD with an average efficiency of 99% and a 
hadron contamination below 1% [1]. 




Fig. 12. HERMES RICH event display for an 
event with a 14.6 GeV electron (right) and 
a 1.5 GeV tt~ (left) in the lower half and a 
5.5 GeV K + in the upper half. See text for 
detailed description. 

lower half, and a 5.5 GeV kaon in the upper 
half of the detector. The solid black points 
mark the PMT hits, while the markers in 
the ring centers indicate the virtual track 
hit points, i.e. the points where the particle 
tracks would intersect the photon detec- 
tor if they were imaged by the mirror. The 
solid lines are spline fits to a few simulated 
photon hits. They indicate where, based on 
the track parameters and the particle type, 
hits could be expected for this event. The 
electron track is easily identifiable as the 
only one with a gas ring; a comparison to 
figure 2 shows that only for the electron a 
gas ring is expected. The momentum of the 
pion is below the pion gas threshold and 
thus it only exhibits an aerogel ring. How- 
ever, the particle clearly must be a pion, 
because the particle momentum is below 
the aerogel threshold for kaons. The kaon 
in the top detector has a well defined aero- 
gel ring, but no gas ring - as is expected for 
a kaon at 5.5 GeV. Any pion at this mo- 
mentum would certainly have produced a 
gas ring, while the aerogel ring for a proton 
would have a much smaller radius. 
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The task of the pattern recognition algo- 
rithm is in principle to associate the vari- 
ous hits in the photon detector of the RICH 
with a ring assigned to a certain track. 
In practice this detailed assignment is not 
needed and only the information relevant 
to the track identification is extracted. The 
relatively low track multiplicity in the HER- 
MES experiment simplifies this task con- 
siderably. 



4-2 Inverse Ray Tracing 

The analysis of the hit patterns is intrin- 
sically complex since the non-linearities of 
the imaging system distort the simple ring 
structure of the emitted light. The influ- 
ence of the imaging system can be removed 
by inverse ray tracing (IRT) [8,9]. In this 
method the Cerenkov angle corresponding 
to each PMT hit is reconstructed from the 
track parameters and the position of the 
PMT. 

The inverse ray tracing problem can be 
formulated as follows: given a track and 
a hit in the RICH photon detector plane, 
at which angle was the photon emitted? 
Assume that the emission point can be es- 
timated. This assumption will be discussed 
below. 

The geometrical problem can be formu- 
lated as follows using the terminology of 
figure 13. Given point E, the likely emis- 
sion point, point D, the detection point 
and C, the center of the spherical mirror 
the photon scatters from, find the point 
S on the surface of the mirror where the 
photon scattered. The properties of point 
S are, in vector notation: 

(1) CS is coplanar with CE and CD. 



D 




Fig. 13. The problem of inverse ray tracing. 

(2) |C75| = R 

(3) <(ES, CS) = <(CS, DS) 

For the mathematical formulation of the 
problem it is easier to switch to an Eu- 
clidean base with C as the origin. The u 

— * 

axis is defined along CE. The v axis is 

— * — * 

coplanar with CE and CD, and oriented 
such that v-CD > 0. The usual caret nota- 
tion is used to indicate a unit vector. The 
third axis (w) is defined such that the base 
defines a right-handed coordinate system. 
The components of the vectors are then: 

CE = (a, 0,0) 
CD = (dcosa, dsino;,0) 
CS = (R cos (3, R sin (3,0) (3) 

Clearly \CE\ = a and \CD\ = d. The an- 
gles are defined in figure 14. 

The coplanarity requirement has been sat- 
isfied by the choice of the axes. The para- 
meterization of S ensures that it will be 
located on the mirror surface. The remain- 
ing parameter (3 now must be determined 
from the equality of the incident and re- 
flected angles. This leads to two (redun- 
dant) equations for the unit vectors: 
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Fig. 14. Definition of the base and the angles 
SC ■ SE = SD ■ SC (4) 

(5) 



SC x SE =SD x SC 



The second equation is a vector equation, 
but since all three vectors are coplanar, the 
only parameter is the length of the vector 
product. 

The first equation reduces to: 

\SD\ d • (cos ol cos (3 + sin o. sin /3) — R 



\SE\ 



a cos (3 — R 



(6) 



The vector equation represents the equal- 
ity of the sines of the angles, and has only a 
non-zero component along the w-axis. This 
component reduces to: 



\SD\ d{— cos a sin/? + sin a cos f3) 



\SE\ asin/3 
Combining the two equations yields: 



(7) 



ad ■ sin(a — 2(3) + 

R- (asin/3 - dsin(a - /?)) =0 (8) 

The solution to the equation can now be 
found from Newton-Raphson iterations. 
Since the result will be the angle (3, no fur- 
ther work is necessary. Using half of the 
angle a as the starting point, convergence 
is usually achieved within three iterations. 



The reflection point on the mirror can then 
be found as follows: 

S = C + ( Rcos P - Rsin l 3 " cosa ) . (f E 
a a ■ sin a 

Rsmd -> , , 
+ , . - CD 9 
asm a 



The inverse ray tracing method reconstructs 
the emission angle for each hit, provided 
the emission vertex is known. For each track 
it is assumed that the hit could be com- 
ing from the aerogel or the gas, and the 
emission vertex is estimated accordingly. 
For each radiator hypothesis, the emission 
angle is then reconstructed. This results 
in a spectrum of reconstructed angles, as 
shown in figure 15, which is clearly depen- 
dent on the radiator hypothesis. Due to 
the very different index of refraction, it is 
rarely a problem to distinguish which radi- 
ator hypothesis is correct. The remaining 
uncertainty in the emission vertex due to 
the finite length of the radiator is unavoid- 
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Fig. 15. Reconstructed angle spectra for the 
assumption that the photon was emitted in 
the aerogel (top) or in the gas (bottom). The 
data are from single track, low background 
electron events with p > 5 GeV. 
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able and results in the contribution A6 em 
to the resolution in table 3. 

Another difference in the treatment of re- 
constructed aerogel and gas angles lies in 
the refraction of aerogel photons at the 
aerogel-gas boundaries. This is approxima- 
tely corrected for under the assumption of 
a flat aerogel-gas boundary and coplanar 
refraction. 



4-3 Likelihood Analysis 



The selection of a most probable particle 
type is made by calculation of the likeli- 
hood that each particle type would gener- 
ate the measured IRT spectrum. Due to 
aerogel edge effects and the finite mirror 
size, the yields of aerogel and gas hits vary 
strongly over the spectrometer acceptance. 
As a result the particle identification is 
largely based on the reconstructed aver- 
age angles, which are less susceptible to 
these effects. The theoretically expected 
angles Q t h are calculated for each parti- 
cle type hypothesis % from the track mo- 
menta. For each particle type hypothesis 
a window is imposed on the spectrum of 
reconstructed angles around this expected 
value. The width of this window (a w ) is set 
to 4 times the single photon resolution o~g. 
However, it is in principle a free parame- 
ter that can be adjusted to optimize the 
performance of the particle identification 
algorithm. The average angle is then de- 
termined from the reconstructed angles in 
the window [9 th - a w ; 9 th + a w \. 

The distribution of measured average an- 
gles for a given particle type and momen- 
tum can be normalized to form a condi- 
tional probability. For known relative par- 
ticle fluxes a particle type probability can 



then be calculated. If the relative parti- 
cle fluxes are not known, the conditional 
probability can itself be used likeli- 
hood. Furthermore, if the average angle 
resolution is assumed to be independent 
of the particle type, the overall normaliza- 
tion of the conditional probability distri- 
bution does not matter. It may be chosen 
so that the maximum value of the likeli- 
hood is 1. With a Gaussian shape of the 
average angle distribution the likelihood is 
then calculated as 



L({9)) = exp 



(10) 



where a^) is the average angle resolution, 
which is calculated from the single photon 
resolution as 



a (9) = 



0-8 



N 



(11) 



The likelihoods for the two radiators, aero- 
gel and g combined to an overall 
likelihood by multiplication. For this pa- 
per, the particle is assigned the type with 
the highest likelihood. Obviously, this re- 
quirement will allow more contamination 
than is acceptable for most physics analy- 
ses. A detailed discussion of the "quality" 
of the particle identification and the result- 
ing tradeoff between contamination and ef- 
ficiency will be the subject of a subsequent 
paper. 



5 Mirror Alignment 



The alignment of the mirror array within 
the detector volume was determined from 
the data. It consisted of three steps: global 
alignment, individual mirror alignment and 
local alignment. 
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Fig. 16. Focusing quality (0}\ for pion gas rings (4.5-5.5 GeV) before and after the individual 
and local mirror alignment. The lines indicate the boundaries of the individual mirrors. 



For global mirror alignment the mirror ar- 
ray was treated as one unit. A three-dimen- 
sional scan of the assumed position of the 
mirror center was performed (in software) 
and the influence on the reconstructed aero- 
gel and gas rings was studied using a set of 
high energy single electron tracks. The per- 
formance was characterized using 24 pa- 
rameters, including the average angles, their 
spread, the single photon resolution, the 
yields and the azimuthal uniformity of the 
gas rings. A careful study resulted in a 
small range of acceptable values for the 
position of the mirror center and a shal- 



low optimum was selected. 

However, both mirror arrays consist of eight 
individual mirrors each, whose alignment 
relative to one another is not necessarily 
perfect. In addition, the construction of 
the mirror allows for local variations of the 
radius of curvature of the surface, espe- 
cially around the stress points introduced 
to compensate for the flexibility of the in- 
dividual mirror segments (positioned at the 
edges of the mirrors). These variations in 
the mirror radius introduce local variations 
in the reconstructed angle. This translates 
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into a larger spread on the reconstructed 
angle, and (if not corrected) a different av- 
erage angle for some parts of the mirror. In 
certain momentum ranges this can lead to 
a misidentification of the particle. 

For this reason in the next two steps of 
the alignment determination the individ- 
ual mirrors are aligned relative to one an- 
other and then each mirror is further sub- 
divided and the pieces are individually ad- 
justed. Steps two and three are based on 
two geometric properties of the mirrors: 
the position of a spherical mirror with 
known radius is unambiguously determined 
by the position of the mirror center. Fur- 
thermore, if it is known how a single pho- 
ton is reflected by the mirror, the mirror 
center can be determined for a given ra- 
dius. 

While the actual track of a photon in the 
detector cannot be determined, it is pos- 
sible to determine approximately how a 
charged particle track would reflect in the 
mirror. This is done by measuring the cen- 
ter of the gas ring created by the track. 
This center can be related to the point 
where a reflected track would hit the pho- 
ton detector by comparing it to an MC 
simulation. For a single track the recon- 
struction of the individual mirror centers is 
still underdetermined. However, it is pos- 
sible to use a large number of tracks with 
gas rings that were reflected on the same 
individual mirror and then calculate the 
optimal mirror center to account for their 
behavior. This optimum is not unique and 
depends somewhat on the starting point. 

For the third step the individual mirrors 
were artificially divided into 3x3 parts 
that were treated as if they were individual 
mirrors. The limiting factor here is the size 
of the gas rings. This 'local alignment' al- 



lows for a better approximation of the real 
mirror surfaces. The result is illustrated for 
the top mirror array in figure 16. The plot- 
ted quantity is the focusing quality 
\6i — (9)\, where 9i is the reconstructed 
Cerenkov angle for a single PMT hit and 
(9) is the average angle for the correspond- 
ing track. This quantity approximately cor- 
responds to a local derivative of the mirror 
surface. Individual mirror and local align- 
ment are based on a large sample of pions 
with momenta between 4.5 and 5.5 GeV. 
They were selected to cover as much as 
possible of the mirror surface while having 
well defined gas rings. 



6 Detector Performance 

The detector performance ultimately is mea- 
sured by how well the various particle types 
are identified. For a given particle identifi- 
cation algorithm this performance is deter- 
mined by the number of detected gas and 
aerogel photons, as well as the single pho- 
ton resolution. Complications arise from 
overlapping rings and background hits. 

6. 1 Yields 

The number of aerogel photons detected 
strongly depends on how much the track 
is affected by acceptance effects from the 
finite size of the mirror array and the tile 
structure of the aerogel radiator. Figure 17 
illustrates the influence of the tile struc- 
ture. The bottom plot is a projection of 
of a horizontal band of tile in the top plot 
between y=30 cm and 36 cm, showing the 
edges of the aerogel tile stacks. The slightly 
lower average yield in the top left tile arises 
from a lower incident electron flux. With 
the given size of the aerogel tiles about 
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Fig. 17. Dependence of electron aerogel yield 
on track position at the aerogel radiator. 
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Fig. 18. Number of fired aerogel (top) and 
gas PMTs (bottom) versus particle momen- 
tum for 'ideal' pion tracks. 

45% of all tracks are affected by the tile 
structure. It would be desirable to develop 
larger aerogel tiles of the same optical qual- 
ity. For example, in the case of 20x30 cm 2 
tiles only about 20% of the tracks would 
be affected. 



Figure 18 shows the number of aerogel PMT 
hits N(aero) and gas PMT hits N(gas) ver- 
sus particle momentum for 'ideal' pion 
tracks that do not suffer from acceptance 
or overlap effects. The two plots show the 
rise towards the asymptotic values of 10 
aerogel hits and 12 gas hits. As a result 
of the acceptance effects, the asymptotic 
number of aerogel hits for all tracks is re- 
duced to about 8. 



6.2 Background 



An inspection of all events from a single 
run (1861 events) showed that roughly one 
third of all particles that produce clear rings 
in the detector do not have a full recon- 
structed track that is associated with them. 
These tracks are low momentum particles 
that do not pass through the entire spec- 
trometer and particles produced in flight 
within the spectrometer. Examples of the 
latter are delta electrons and electron-posi- 
tron pairs. These 'trackless rings' act as 
background for the rings that are associ- 
ated with tracked particles. In events with 
more than one track in one detector half 
the rings of the tracked particles of course 
are 'background' for one another in the 
same way. 

In addition to the trackless rings, there are 
several sources of background photons that 
do not necessarily result in ring structures. 
These include Rayleigh scattered photons, 
Cerenkov photons produced in the lucite 
window, proton beam correlated back- 
ground showers that hit the PMT matrix 
directly and scintillation in the gas. The 
electronic and PMT noise in the detector 
is a very small effect, amounting to only 
about 1 fired PMT every 5 events. 
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6.3 Average Angles and Resolution 



As the likelihood analysis is based on the 
average angles, it is of particular impor- 
tance that the angles are correctly recon- 
structed and that their momentum depen- 
dence is understood. Figure 19 shows the 
reconstructed average aerogel angle for pi- 
ons. The data were fit with a theoretical 
curve with the aerogel index of refraction 
as the only free parameter. The resulting 
curve and index of refraction (n= 1.0304) 
are in excellent agreement with the theo- 
retical expectation and the optically mea- 
sured index of refraction. Figure 20 shows 
the reconstructed gas angles for pions to- 
gether with the theoretical curve for an in- 
dex of refraction of 1.00137. The system- 
atic difference between data and theoret- 
ical curve at low momenta is due to the 
finite size of the PMTs. 

The MC simulation of the RICH detec- 
tor plays an important role in the under- 
standing of efficiency and contamination 
values for the RICH particle identification. 
Hence, it is necessary that the detector re- 
sponse is adequately described by the sim- 
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ulation. Figures 21 and 22 show the com- 
parison of reconstructed aerogel and gas 
angles for ultra-relativistic electrons 
(p > 5 GeV). 

For the comparison the data and the MC 
results were normalized to the same max- 
imum value. 

Figure 21 shows the excellent agreement 
in the central part of the distribution. The 
wider tails of the experimental distribution 
are due to several effects that will be dis- 
cussed below. The flat background past the 
tails is another important difference be- 
tween the MC and the data. The differ- 
ence in the aerogel spectra is due to the 
presence of additional contributions to the 
resolution, which are listed in table 4. A 
significant effect comes from the aerogel 
tiles, which have sharply curved surfaces 
near the edge. This refracts the photon in 
an unpredictable way, which gives rise to 
the A9 ti i e contribution. It also contributes 
to the background. A9 tile was estimated 




Fig. 19. Reconstructed average aerogel an- 
gle versus particle momentum for pions. The 
solid line represents a fit with n=1.0304. 
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Fig. 20. Reconstructed average gas angle 
versus particle momentum for pions. The 
solid line represents a theoretical curve with 
n=1.00137. 
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Fig. 21. Normalized distributions of recon- 
structed aerogel angles for single, low back- 
ground electrons (p > 5 GeV). Histogram: 
MC, points: data. 
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Fig. 22. Normalized distributions of recon- 
structed gas angles for single, low background 
electrons {p > 5 GeV). Histogram: MC, 
points: data. 

from comparing the resolution for tracks 
that are not affected by the tile edges to 
those that are. A9 nvar accounts for varia- 
tions in the index of refraction between the 
aerogel tiles. Its value is determined from 
the measured indices of refraction. There 
is also a 'forward scattering' contribution 
A6f w [21], which must be taken into ac- 



count, but which, unlike the contribution 
A6 c hr, is not included in the Monte Carlo 
simulation. 

For the reconstructed gas angles in figure 
22, the width and shape of the distribu- 
tions is basically in good agreement. The 
MC distribution is slightly wider and less 
'Gaussian' than the experimental data. As 
this is the distribution of the individual 
reconstructed angles (with one entry per 
fired PMT), this is likely an effect of the 
MC overestimating the photoelectron yield 
for the gas. Due to the non-linear relation- 
ship between the number of photoelectrons 
and the number of fired PMTs the addi- 
tional photoelectrons add relatively more 
to the sides than to the center of the dis- 
tribution. 

Pressure and temperature fluctuations 
change the density of the gas, and there- 
fore its index of refraction. This contri- 
bution, A8 press , was determined from the 
observed maximum pressure fluctuations. 
This error can be eliminated by taking the 
measured pressures into account in the re- 
construction. 





aerogel 


C4F10 


AOtile 


~3.0 mrad 




AOnvar 


~1.1 mrad 




A9f w 


~0.9 mrad 




A6p T ess 




~1.0 mrad 


A6 e xp 


7.6 mrad 


7.5 mrad 


AQ m irr 


~2.3 mrad 


~2.2 mrad 



Table 4 

Additional contributions to the single photon 
resolution in aerogel and C4F10. The contri- 
bution to the resolution related to the mirrors 
(A9 m i rr ) has been obtained by subtracting 
the other relevant contributions from A9 exr> . 
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The experimental value of the total an- 
gular resolution for single photons, A8 exp , 
was determined for high energy single lep- 
ton tracks that are not affected by any 
tile edge or mirror acceptance effects (ta- 
ble 4). The data were taken from a period 
of stable atmospheric pressure and events 
with a large background were excluded. 
The reconstructed angle spectra were fit- 
ted with a Gaussian plus a linear back- 
ground. As a result of the selection of these 
'ideal' tracks the aerogel tile effects (A8ta e ) 
do not contribute to A9 exp for aerogel and 
A9 press does not contribute to A9 exp for 
gas. The experimental gas resolution there- 
fore is a combination of the effects that 
are collected in A9mc (cf. table 3) and the 
hitherto unknown contribution of the mir- 
ror imperfections A9 mirr . This contribu- 
tion in particular refers to the non-sphericity 
of the mirror array and the diffuse reflec- 
tion component of its surface. Under the 
assumption that the quadratic sum of A9mc 
and A9 mirr yields the experimental value 
A9 exp for gas, the mirror contribution is 
calculated to be about 2.2 mrad. In the 
case of the aerogel also A9 nvar and A9f w 
must be included. The assumption that A9 ex 
for aerogel is the quadratic sum of these 
contributions with A9 M c (cf- table 3) leads 
to a value of 2.3 mrad for A9 mirr , in good 
agreement with the value extracted from 
the gas angles. 

A comparison of tables 3 and 4 shows that 
the dominant contribution to the resolu- 
tion comes from the pixel size (A9 pix ). 



6.4 Efficiencies 

Monte Carlo predictions for the efficien- 
cies and purities of the hadron identifica- 
tion using the IRT likelihood algorithm are 



given in table 5. The MC simulation of 
the RICH detector used for this predic- 
tion contains a detailed simulation of the 
aerogel radiator geometry as well as exper- 
imental background data extracted from 
trackless events. The data are from 1.84 
million MC events using a deep inelastic 
scattering generator. The particles are iden- 
tified as the particle type with the high- 
est IRT likelihood. As mentioned before, 
stronger constraints on the particle likeli- 
hoods will improve the purity of the sam- 
ple, albeit at the expense of the efficiency 



id. as 


7T 


K 


P 


purity 


7T 


619390 


8050 


23008 


0.95 


K 


17282 
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Tabic 5 



Typical efficiencies and purities for the RICH 
hadron identification based on the IRT likeli- 
hood analysis of 1.84 million DIS MC events, 
p > 2 GeV. 



Tab: 



Decay 


Efficiency 


p° — ► 7T + 7T~ 


e n = 0.915 ±0.024 


K s — ► 7T + 7T~ 


e n = 0.900 ± 0.005 


(j) -► K+K- 


e K = 0.750 ±0.007 


A — ► pir~ 


e p = 0.726 ± 0.010 




e n = 0.890 ±0.011 


le 6 



Typical momentum integrated efficiencies de- 
termined from decaying particles based on the 
IRT likelihood analysis. 

It is not possible to use another detector to 
create clean hadron samples to study the 
detector performance. However, it is possi- 
ble to use decaying particles for the same 
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purpose. Samples of p, and K s mesons 
as well as A hyperons were used to deter- 
mine the identification efficiencies for pi- 
ons, kaons and protons. The momentum 
integrated results (p>2 GeV) are shown in 
table 6. The MC values in table 5 and the 
experimental values in table 6 cannot be 
compared directly, because the efficiencies 
depend strongly on the particle momen- 
tum as well as on the event topology. 



p S 1 
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0.92 
P * K 0.06 
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Fig. 23. Identification probabilities that a 
pion is identified as pion, kaon, proton or not 
identified (X) for single tracks per detector 
half. DIS MC simulation (circles) in compar- 
ison with experimental data from p°-decays 
(solid squares); both results are based on the 
IRT likelihood analysis. 

This dependence on the event topology can 
be excluded by comparing subsamples with 
only one track per detector half. Figure 23 
compares MC and experimental data for 
the case of single track pions. The four 



plots show the identification probabilities 
that a pion is identified as pion, kaon, 
proton or not identified (X). The pion iden- 
tification efficiency P% for single tracks is 
above 92% even at low momenta, and above 
98% for momenta larger than 4 GeV. The 
MC data are taken from the sample also 
used for table 5, while the experimental 
data are from p°-decays. MC and experi- 
mental results agree well,with the possible 
exception of the lowest moment bin. 

The comparison of MC and experimental 
data becomes more involved for event topolo- 
gies with overlapping rings and as the open- 
ing angle for 0-decays in HERMES is small 
there are basically no experimental data 
with one single kaon per detector half. A 
detailed evaluation of the detector perfor- 
mance that takes the event topology and 
momentum dependence into account in de- 
tail will be included in a forthcoming pa- 
per. 



7 Conclusion and Outlook 



The HERMES RICH detector has been 
constructed and installed within 13 months 
after its approval. It has been operating 
routinely as part of the HERMES experi- 
ment since its installation in May 1998. Its 
operation has been stable and reliable for 
more than two years. The single photon 
resolution for ideal tracks is close to the 
MC predictions. The particle identification 
based on the inverse ray tracing technique 
has been implemented and its likelihood 
analysis has been optimized. More elabo- 
rate particle identification schemes are un- 
der development. Hadron identification by 
the RICH detector will be a crucial feature 
of the analysis of the current (1998 - 2000) 
and future HERMES data. 
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